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Abstract 

Local kinetic and chemical equilibration is studied for Au+Au collisions at 10.7 
AGeV in the microscopic Ultrarelativistic Quantum Molecular Dynamics model 
(UrQMD). The UrQMD model exhibits dramatic deviations from equilibrium dur- 
ing the high density phase of the collision. Thermal and chemical equilibration of 
the hadronic matter seems to be established in the later stages during a quasi- 
isentropic expansion, observed in the central reaction cell with volume 125 fm 3 . For 
t >10 fm/c the hadron energy spectra in the cell are nicely reproduced by Boltz- 
mann distributions with a common rapidly dropping temperature. Hadron yields 
change drastically and at the late expansion stage follow closely those of an ideal gas 
statistical model. The equation of state seems to be simple at late times: P = 0.12e. 
The time evolution of other thermodynamical variables in the cell is also presented. 

Key words: Statistical model, equation of state, thermalization, Monte-Carlo 
model for relativistic heavy ion collisions. 
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1 Introduction 



The main goal of the relativistic heavy-ion experiments at Brookhaven and 
CERN is to study the properties, e.g. the equation of state (EoS), of strongly 
interacting hot and dense hadronic matter produced in the course of nuclear 
collisions. A question of great importance in these studies is to determine the 
EoS and whether local (LTE) or even global thermodynamical equilibrium is 
reached in the system or not. This is a crucial point for a large group of models, 
widely used for data analysis. Despite the long history of such investigations 
(e.g. [1-7]), this question remains still open. 

One can try to check the assumption of LTE by analyzing the experimental 
particle number ratios and hadron momentum spectra (e.g. [8-11]). Due to 
the time duration of chemical (inelastic) and kinetic (elastic) freeze-out of the 
system in A+A collisions [12-15], the answer is not conclusive. In contrast to 
the microscopic approaches [16-20], macroscopic models like hydrodynamical 
[21-23] or thermal [7,9,10] models adopt LTE as ad hoc assumption, use the 
Equation of State [24] as a parametrization and apply usually an instanta- 
neous freeze-out of hadrons, produced from many different space-time cells. 
The freeze-out conditions for these cells - temperatures, baryonic chemical 
potentials and collective velocities - are actually quite different and model de- 
pendent (see e.g. [11,15,25]). Also note that the experimental particle spectra 
in elementary e + e~ and hadron-hadron collisions at high energies already look 
"thermalized" [25,26]. 

Because of these difficulties it is very promising to check the LTE and to 
extract the EoS at different stages of relativistic heavy ion reactions based 
on a detailed microscopic analysis. Here, we employ for a such analysis the 
microscopic Ultrarelativistic Quantum Molecular Dynamics model (UrQMD) 
[16]. 



2 Ultrarelativistic Quantum Molecular Dynamics 

UrQMD is a N-body transport model designed for describing heavy ion colli- 
sions in the laboratory energy range from several MeV to several TeV per nu- 
cleon. The detailed presentation of UrQMD and underlying concepts together 
with a comparison with experimental data is now publically accessible [16]. 
The model treats binary elastic and inelastic hadronic collisions and many- 
body resonance decays. The inelastic collisions and decays are the only source 
for changing the chemical composition of the system, while elastic collisions 
change only the momentum distributions of the hadrons. Both processes drive 
the system towards thermodynamical equilibrium. However, the time scales 
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may be too short in heavy ion collisions for actually achieving LTE. 

The UrQMD model is based on string and resonance excitation both in pri- 
mary nucleon-nucleon collisions from projectile and target, some also in sec- 
ondary hadronic collisions. There are 55 baryon and 32 meson states as dis- 
crete degrees of freedom in the model as well as their anti-particles and explicit 
isospin-projected states with masses up to 2.5 GeV. From 1.5 GeV the strings 
can be populated as continuous degrees of freedom. The experimental hadron 
cross sections and resonance decay widths are used when possible. At high 
energies a string mechanism is implemented to simulate soft hadronic interac- 
tions. Hadrons, produced through the string decays, have non-zero formation 
time which depends on the four-momentum of the particle. Newly produced 
particles cannot interact during their formation time. Leading hadrons (con- 
taining the constituent quarks) interact within their formation time, but with 
reduced cross sections, proportional to the number of original valence quarks. 
The Pauli principle is applied to BB collisions by blocking the final state if the 
outgoing phase space is occupied. No Bose effects for mesons are implemented 
in UrQMD. 



3 UrQMD analysis of Au+Au collisions. 

Now let us use the UrQMD model to explore to which extent LTE can be 
reached at least in a small cell around the origin of central (b=0 fm) Au+Au 
collisions at 10.7 AGeV. The formation of hot and dense hadronic matter is 
expected in this reaction, therefore the comparison of local system properties 
predicted by UrQMD with those of the ideal hadron gas looks promising. 

For our analysis we consider a cubic cell of volume V = 5x5x5 fm 3 centered 
around the origin of center-of-mass of Au+Au system. The collective velocity 
of this cell is equal to zero because the system geometry gives no preferable 
direction for the collective motion. The cell size should be small enough to 
avoid collective flow (and streaming) and not to underestimate densities inside 
the cell. On the other hand it should be large enough to contain enough 
particles inside the cell, which guarantees reasonably small fluctuations of 
particle observables in the cell. From our analysis (see below) it appears that 
a cubic cell with a side of 5 fm satisfies simultaneously both above requirements 
for 200 Au+Au central collisions at AGS energy. 

We start by analyzing the time evolution of different physical quantities in the 
Au+Au center-of-mass frame. Time t — fm/c corresponds to the moment 
when the two Lorentz contracted nuclei touch each other. The maximum over- 
lap of the nuclei occurs at about 2.5 fm/c. At time t = 2Ra u / (jcmVcm) —5.5 
fm/c the freely streaming nuclei would have passed through each other. 
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The reaction is dominated by strong momentum anisotropies for all particles 
inside the cell. During the entire high energy density phase, £ <10 fm/c, the 
widths of the velocity distributions in the longitudinal (z—) direction, a z , are 
much larger then the widths in the transverse (x—, y—) directions, a x and 
<j y , for each hadron species. However, for £ >10 fm/c the anisotropies inside 
the central cell disappear and the widths oi become approximately equal. The 
velocity component distributions of nucleons in the cell, shown in Fig. 1 at two 
different times, confirm this result. The shapes of the velocity distributions 
are very different for v z and v± at t—5 fm/c. They become very similar at 
£=10 fm/c. This approximate isotropy of hadronic velocity distributions in the 
central cell appears to be due to the escape of fast or non-scattered particles 
from the cell and as a result of many hadronic rescatterings. At £ =10 fm/c 
the newly produced particles in the cell have undergone from 1.5 (pions and 
K) to 5 (K) elastic collisions, while the baryons have suffered more than 20 
strong interactions. Fig. 2 shows the large difference between the hadrons 
containing original constituent quarks (< n co u > = 10 — 30) and newly bred 
mesons (< n co u >= 1 — 3). There are still some collisions in the cell at £ >18 
fm/c. Fig. 3 shows that for £ <18 fm/c only a small fraction of the particles 
in the cell is frozen out before that time. 

The isotropy of the momentum distributions for all hadron species in the 
cell is a necessary prerequisite for LTE. The next observation is that hadron 
energy spectra in the central cell at £ > 10 fm/c are nicely fitted by Boltzmann 
distributions, exp(— Ei/T), with rapidly falling T— values vs. time. At each 
moment, however, all particle species % exhibit practically the same slope. The 
"temperatures", T, agree to within AT =10 MeV. 



4 Ideal Gas Statistical Model. 

To study chemical equilibrium in detail, let us now compare the hadron spectra 
and yields in the cell, as predicted by UrQMD model, with those of a Statistical 
Model of an ideal hadron gas (SM) [27,28]. For this comparison of the UrQMD 
results to the SM model the time interval from 10 to 18 fm/c is chosen. The 
largest values of the energy density, e = 21e , and baryonic density, pb = 10po, 
in the central reaction zone are reached earlier, at £ =4-5 fm/c. At this time, 
however, the system of particles in the cell is not yet thermalized and still 
exhibits strong velocity anisotropy. At £ = 4 — 5 fm/c there is also a large 
fraction of non-formed and non-interacting particles in the cell, particularly 
mesons (see Fig. 3). Hence, the physical interpretation of this early stage of 
the reaction should - in the present approach - not be carried forward in a 
SM model. For £ >10 fm/c the fraction of non-formed particles goes to zero 
rapidly. On the other hand the total number of particles in the cell is small at 
£ > 18 fm/c and the system approaches its thermal freeze-out (no collisions 
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occur) stage. In this paper we restrict our analysis of LTE to the time interval 
between t=10— 18 fm/c. We do not discuss the "initial" stage with high energy 
density and the "final" kinetic stage of particle freeze-out. 

For a thermalized ideal gas the momentum spectrum of particle species "i" in 
the rest frame of the cell can be presented as 

dNi Sl V fnA ( Ei\ 9iV 

where V is the volume of the cell, T is the temperature of the system in the 
cell, g iy TOj, Ei = (p 2 + m 2 ) 1 / 2 and pi are the degeneracy factor, the mass, the 
energy and the chemical potential of the hadron species "i" , correspondingly. 
Here, Boltzmann distribution functions, /j, are used for all hadrons. Quantum 
statistical effects are not included in the present analysis. 

With the term thermodynamical equilibrium we specifically mean both kinetic 
and chemical equilibrium. The chemical potential of particle i in LTE can be 
represented as 

IH = + Si/is (2) 



in terms of the baryonic-, /ib, and strange-, chemical potentials {pi = 0, ±1, 
Si = 0, =Fl, =f2, =1=3 are, respectively, the baryonic- and strange- charges of the 
particle species "i"). The electric chemical potential considered in [29,30] is 
neglected. 



5 Results and discussion. 



Here, 200 Au+Au events, generated with UrQMD, are analyzed using the 
following procedure. The average energy density, e, baryonic density, ps, and 
net strangeness density, ps, in the cell are calculated for 10< t <18 fm/c. The 
contributions of both formed and not yet formed particles are included. The 
values of e, ps and p s are inserted in the l.h.s. of the following equations: 

£ = J2 £ i ' PB = h i U i ' PS = J2 S i U i ■ ( 3 ) 

i i i 



Here 
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are the number- and energy density of particle species "i" predicted by the 
ideal gas model. The ideal gas model employs the same hadron species "i" as 
the UrQMD model. Solving the system of three equations (3) for the three 
unknown parameters, T, p B and p,s, we find the results presented in Table 1. 



5.1 Thermalization in central reaction zone. 

Now we can turn to the predictions of the Statistical model (SM) for the 
specific single particle spectra, Eq. (1), by using the intensive parameters T, p# 
and fis, extracted from the average energy density, e, baryon density, p#, and 
strangeness density, ps- Also the pressure and the entropy can be calculated 
swiftly (see below). The Boltzmann energy spectra Eq. (1) predicted by the 
SM with T, p B , Ps values from Table 1 and the spectra obtained in the UrQMD 
model in the central cell for different particle species agree surprisingly well 
for all instants t between 10 and 18 fm/c. Fig. 4 shows this comparison for 
t=13 fm/c. Between 10 and 18 fm/c also the hadronic yields, obtained in the 
SM, are rather close to the UrQMD data, see Fig. 3. These results seem to 
indicate that the LTE is closely resembled for the hadronic matter in central 
cell of central Au+Au collisions at 10.7 AGeV at least for t =10 to 18 fm/c. 
The statistical model parameters, i.e. the main thermodynamic characteristics 
of the cell change rapidly with time. This clearly demonstrates that a fireball 
type description of hadronic matter is inadequate. 

Table 1 shows that during this time interval (10 - 18 fm/c) the baryonic density 
and temperature in the cell change from p B = 2po to 0.3po (po — 0.16 fm -3 
is the normal nuclear density) and from T =150 MeV to 100 MeV. The time 
evolution of the central cell in (T, p#— ) plane is shown in Fig. 5. It resembles 
closely an isentropic hydrodynamical expansion. Let us develop this idea in 
detail. 



5.2 Equation of State. 

The hadron pressure is given in the statistical model by 

oo 2 

Using Eq. (4) we find that the equation of state, P(e), of the hadronic matter 
in the central region of Au+Au collisions at 10.7 AGeV is nearly linear with 
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e. In fact, P(e)/e is constant for the whole time interval £=10-18 fm/c 



P/e = 0.12 . (5) 



This constant corresponds to the square of the speed of sound in the medium, 
c 2 s = (dP/de) =0.12, which is very similar to that in a resonance gas, c 2 s =0.14, 
obtained in [4]. The constant behaviour of P(e)/e confirms that transport 
models like UrQMD and RQMD do not predict any phase transition which 
would be manifest in the softest point of the EoS, P(e)/e [11]. 



In Fig. 6 we compare the ideal gas pressure Eq. (4) with longitudinal, P{™ c 
and transverse, P{™^}, pressures obtained in UrQMD as [31,32]: 



3V(ml + pl)W • 



pmic \p Ph{x,y,z} 

r {x,y,z} ot//™2 i „2M/9 ' V U / 



where V is volume of the cell, ph represents the particle momentum and the 
sum in Eq. (6) is taken over all hadrons h in the cell. Fig. 6 shows a strong 
difference between P{™ c and P^ l y} f° r £ <10 fm/c. That this is due to the 
longitudinal flow in the cell at early times can be seen from the dashed curve, 
obtained for a cell with a shortened longitudinal side, Az < 1 fm. The differ- 
ence between P± and Py vanishes already at £ =6 fm/c, but at so early times 
there is still no equilibrium in the cell: the slopes of the energy spectra are 
quite different for different particle species. Only at £ =10 fm/c the hadronic 
matter in the cell becomes equilibrated. At this time the microscopic pressure 
Eq. (6) in both cells also becomes nearly isotropic and approximately equal to 
the ideal gas one Eq.(4), 3Pg c = 3P { ™J } ^ p. The strong anisotropy of the 
transverse and the longitudinal components of the pressure is also observed 
in three-fluid hydrodynamical model for the initial stages of central Au+Au 
collisions at AGS [33]. This finding is in stark contrast to the one-fluid hydro- 
dynamical model, where the pressure is assumed to be locally isotropic at all 
times. 

Corse graining problems inhibit the direct microscopic calculations of the en- 
tropy. However, the entropy density s can be calculated from the ideal gas 
model (using the microscopic e, ps, Ps as input, see above) by the thermody- 
namical identity 

e = Ts + pbPb + PsPs - P- 

The strange particle contribution is negligibly small (psPs ~ 1 MeV maxi- 
mally) in comparison to the other terms which are of order 100 MeV. The 
entropy per baryon in the cell is nearly constant and is equal to 

S/ A = s/p B 12 . (7) 
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Hence, the time evolution of the cell is nearly isentropic for times £=10-18 fm/ c 
and, therefore, similar to an ideal hydrodynamic expansion. 



5.3 Comparison with infinite equilibrated matter. 

Quantum statistical effects in the ideal hadron gas are small for the temper- 
atures and chemical potentials, T and /i, from the Table 1. Bose statistics 
gives the strongest effect for the pions - the Bose pion number density at T— 
100-150 MeV is about 10% larger than the Boltzmann approximation. The 
quantum statistical effects are only 1-2% for other hadrons. 

To make a conclusion about LTE in heavy ion collisions in a given micro- 
scopic model one has to perform a comparison not only to the Statistical 
Model, but also to an equilibrated infinite matter system within the same mi- 
croscopic model. The energy spectra and the particle yields in the central cell 
of Au+Au collisions were compared to the equilibrated hadronic matter cal- 
culations obtained in a box with periodic boundary conditions by employing 
the UrQMD model to such a closed system with the approximate microcanon- 
ical infinite matter. This closed box calculations must clearly be distinguished 
from the open system, off equilibrium, in the central cell of heavy ion col- 
lisions discussed before. The particle yields and energy spectra predicted in 
the UrQMD box calculations [34] for each set of e, pb and p$ values taken 
from Table 1 for the different time steps are in a good agreement with the 
UrQMD predictions in the central cell of Au+Au collisions (see Fig. 4). This 
additionally confirms that at the intermediate stage of the reaction t =10-18 
fm/c the hadronic matter in the cell is in thermodynamical equilibrium. This 
also means that the closed box calculations (as well as SM) can be used to 
describe - of course not to predict - the time evolution of locally equilibrated 
open system of the heavy ion collisions discussed above. 

The static box values agree approximately with the ideal gas hadron model. 
However, there are enhancements of the low-energy pion yields as compared 
to the ideal gas results observed in both the UrQMD cell- and the box- calcu- 
lations, seen in Fig. 4. The box calculations [34], at ps = Po, show that this 
enhancement becomes even stronger with increasing energy density. This en- 
hancement is due to the violation of the detailed balance conditions for string- 
and many-body resonance decays in microscopic models like UrQMD [34] . The 
additional pion production from many-body decay channels acts like a posi- 
tive pion chemical potential. The Bose enhancement of pions should therefore 
be more pronounced at higher energy densities [35]. 
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6 Conclusions. 



Equilibration of hadronic matter produced in the central cell of Au+Au colli- 
sions at 10.7 AGeV has been studied in the microscopic UrQMD model. It is 
found that 

• There are large deviations from thermal equilibrium for the high density 
phase, t <10 fm/c. 

• Local thermal and chemical equilibration seems to set in, in the central cell, 
for the expansion stage, £=10-18 fm/c. 

• Both the different particle yields and the energy spectra in the central cell 
change drastically with time. Still, their instantaneous values agree, sepa- 
rately, for different time steps t >10 fm/c, with correspondingly adjusted 
ideal hadron gas model values. 

• The temperatures and baryonic densities extracted in the central cell drop 
rapidly from T =150 to 100 MeV and from p B = 2p to 0.3p for t =10 to 
18 fm/c. 

• The entropy per baryon is nearly constant, S/A = s/p B = 12, during the 
time interval t =10-18 fm/c. 

• The equation of state of this complex mixture of different particle "fluids" 
has a very simple form P(s) = 0.12e for this range of T, p B , ps~ values. 

• A substantial enhancement of low-energy pions is observed in the UrQMD 
calculations as compared to the ideal gas results. This effect increases at 
higher energy densities due to the string- and many-body decays of the 
resonances. 

• At t= 10-18 fm/c the spectra and hadron multiplicities in the central reac- 
tion cell are in agreement with those found for equilibrated infinite matter in 
UrQMD calculations in a box with periodic boundary conditions. The time 
scale of the thermal and chemical equilibration is, however, much shorter in 
the present calculations for nuclear reactions in the cell than in the UrQMD 
box calculations. It is the short-lived pre-equilibrium high density initial 
phase which prepares the extreme entropy densities forming doorway state 
to the quasi-isentropic expansion for t >10 fm/c. 

Let us stress again that these studies support a picture of strong deviation from 
local thermodynamical equilibration in the central cell for t <10 fm/c in A+A 
collisions at AGS energies. Large collective streaming-through flow is found 
even in the central cell at early times. The whole system never reaches global 
thermodynamical equilibrium. However, local equilibration is approached in 
central cell for the later times, t >10 fm/c. These flow effects are not large, 
however, at t >10 fm/c in the central cell, that is why it was chosen. The 
averaged squared velocities in the central cell at t=10 fm/c are ({v{ low ) 2 ) = 
0.04 and ((v f J™) 2 } ^ 0.01 for all particle species. These values are much 
smaller than the thermal velocities of the various hadrons in the cell at that 
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time. For nucleons one finds (vf h ) = 3T/rriN = 0.5. Thermal velocities for pions 
and kaons are even larger. Our analysis neglects this small collective motion 
inside the central cell at t= 10-18 fm. More careful studies give corrections to 
this picture: in Fig. 4 one observes that the "best fits" to the particle spectra 
would give Tjv,a > T^k' small collective motion inside the cell can imitate 
slightly larger (several percents) effective temperatures for heavier particles in 
the ideal gas calculations. Much larger flow velocities are found for all cells 
other than the central one. 

A similar analysis has been done at lower energies in [31,36,37]. The next step 
is to study the approach to LTE for the earlier phase of the evolution and 
at higher collision energies (from the AGS to the SPS). This requires much 
smaller cell sizes and correspondingly larger samples of A+A events. 
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Table 1 

The thermodynamic characteristics of the hadronic matter in the central cell 
(V=125 fm 3 ) in central Au+Au collisions at 10.7AGeV. The energy density, e, 
baryonic density, ps, and net strangeness density, ps, are obtained from the micro- 
scopic UrQMD calculations. The temperature, T, baryonic chemical potential, ps, 
and strange chemical potential, ps, are extracted from the ideal hadron gas model, 
using e, ps and ps as input. 



t 


e 


PB 


Ps 


T 


PB 


PS 


(fm) 


(MeV/fm 3 ) 


(fm- 3 ) 


(fm" 3 ) 


(MeV) 


(MeV) 


(MeV) 


10 


675 


0.33 


-0.008 


147 


510 


128 


11 


511 


0.26 


-0.007 


140 


519 


122 


12 


384 


0.20 


-0.006 


134 


527 


116 


13 


293 


0.16 


-0.004 


128 


534 


112 


14 


222 


0.12 


-0.003 


122 


542 


107 


15 


170 


0.10 


-0.003 


116 


553 


101 


16 


134 


0.08 


-0.003 


111 


560 


96 


17 


102 


0.06 


-0.003 


106 


567 


87 


18 


80 


0.05 


-0.002 


101 


574 


79 
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v x , v y , v z 

Fig. 1. Nucleon velocity distributions dN/dvi (i = z (•), x (□) and y (O)) i n 
central cell of Au+Au collisions at 10.7 AGeV at t =5 (upper frame) and 10 (lower 
frame) fm/c. The solid line corresponds to the Maxwell-Boltzmann velocity distri- 
bution dN/dvi ~ exp(— m^vf /2T) for the "non-relativistic" (vi < 0.6) nucleons 
with temperature T = 147 MeV as obtained in ideal gas model. 
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Fig. 2. The number of collisions, < n co u(t) >, in the central cell of Au+Au collisions 
at 10.7 AGeV versus time, t, for nucleons, deltas, lambdas plus sigmas, pions, kaons 
and antikaons. 
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Fig. 3. The number of particles in the central cell of Au+Au collisions at 10.7 
AGeV as a function of time as obtained in UrQMD model (histograms). Solid lines 
correspond to all hadrons in the cell, dashed lines - to interacting particles and 
dotted lines - to formed hadrons. The numbers of frozen particles in the cell are 
shown by dot-dashed lines. The points represent the predictions of the ideal gas 
model as calculated at the same energy, baryonic and strangeness densities as for 
the UrQMD central cell. 
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Fig. 4. Energy spectra of iV (•), A (□), ir (A), K + (y) and A (*) in the central 
125 fm 3 cell of Au+Au collisions at 10.7 A-GeV at t=13 fm/c (black points) and in 
box calculations (open points) are fitted by Boltzmann distributions, Eq. (1) (lines) 
using the parameters T=128 MeV, ^#=534 MeV, = H2 MeV, as obtained in 
the ideal gas model. 
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Fig. 5. The evolution of the baryon density, ps, and temperature, T, in the central 
cell of Au+Au collisions at 10.7 AGeV at times t=10-18 fm/c. The entropy per 
baryon is constant s/ps ~ 12 thus exhibiting quasi-isentropic expansion in this 
central cell. The equation of state has the form P « 0.12e. 
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Fig. 6. The longitudinal (3 • P{ z }, solid histogram) and the transverse (3 • P{ x ^ y ), 
dashed curves) diagonal components of the microscopic pressure tensor in the central 
cell (Ax x Ay x Az = 5x5x5 fm 3 (histograms) and 4x4x1 fm 3 , (dashed line)) 
of Au+Au collisions at 10.7 AGeV calculated from the virial theorem Eq. (6) are 
compared to the ideal gas pressure Eq. (4) (solid line). 
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